Tight coupling in thermal Brownian motors 
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We study analytically a thermal Brownian motor model and calculate exactly the Onsager co- 
efficients. We show how the reciprocity relation holds and that the determinant of the Onsager 
matrix vanishes. Such condition implies that the device is built with tight coupling. This explains 
why Carnot's efficiency can be achieved in the limit of infinitely slow velocities. We also prove 
that the efficiency at maximum power has the maximum possible value, which corresponds to the 
Curzon-Alhborn bound. Finally, we discuss the model acting as a Brownian refrigerator. 
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Can Carnot's efficiency be achieved by a thermal Brow- 
nian motor? This question is of great importance, not 
only for theoretical statistical physicists but for the tech- 
nological construction of micro and nano machines. In 
fact, the practical and relevant issue is to investigate the 
instrinsic features with which a device should be built in 
order to operate optimally at the limitations established 
by first principles. Here we study a theoretical model 
for a Brownian motor and show that it has the so called 
tight coupling property. This gives a fundamental ex- 
planation for the fact that the system has Carnot and 
Curzon-Alhborn's efficiency bounds. 

Let us consider the model proposed by Sakaguchi p] 
for the Brownian motion of a Feynman-like ratchet device 
. It consists on a Langevin equation in the overdamped 
regime that accounts only for one degree of freedom at 
temperature T. The particle moves under the action of 
a spatially periodic and asymmetric potential V(x) and 
an external force F. The equation of motion is 



-V'(x)-F + r ] (t), 



(1) 



where the friction coefficient has been taken equal to one 
by redefining the time scale. The thermal bath is mod- 
eled by a Gaussian white noise rj(t) with the usual fluctu- 
ation dissipation relation < r)(t)rj(t') >= 2T6(t — t'). We 
take ks = 1, which fixes the energy units. The ratchet 
potential is illustrated in Fig. QJ. 

The second thermal bath at temperature To (which 
is actually necessary to break detailed balance) is intro- 
duced through the following boundary condition for the 
steady probability distribution, 



P(l) = P(Q) e - v/To . 



(2) 



This condition ensures the correct thermodinamical 
properties of the model and incorporates the transition 
probabilities suggested by Feynman, namely, that the 
probability to overcome the barrier in the forward di- 
rection (from left to right) is proportional to e -(Y+Fi)/T 








Figure 1: Scheme of the ratchet potential V(x) and temper- 
ature profiles in Sakaguchi's model. 



(from the Langevin equation Q), while the probabil- 
ity to cross the barrier backwards is e~ v ^ T ° (from the 
boundary condition (2J). Sakaguchi's model can also be 
understood as a Brownian particle moving alternately in 
hot and cold reservoirs along space (see Ref. Q) in the 
limit in which the cold region tends to zero. 

The condition J2J implies that P{x) does not connect 
with continuity the gap points, in which V'(x) is infi- 
nite. From Eq. |Q, one can write the equation for the 
probability distribution in the steady state, 



(V'(x) + F)P{x) + TP'{x) = -J, 



(3) 



where J is the constant and uniform probability flux. 
This equation is linear and can be solved analytically 
with two unknown constants, J and P(0), that are eval- 
uated imposing normalization to P(x) and the boundary 
condition J3J) . The observable of interest is the mean ve- 
locity and it is evaluated in terms of the flux, already 
obtained, as 



v = (x) =Jl = 



a 2 T 



e- a - e- b 



I a(e- b - e~ a ) + (1 - e~ b )(l - e~ a ) 
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(4) 

where a = (V + IF)/T and b = V/Tq. This expression is 
exactly the same that the one derived in Ref. 0] if the 
limit L2 = is performed. 
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Figure 2: Energetics of a thermal motor: heat flow Q from 
the hot source at T gets into the system S and mechanical 
power W is performed at the same time that heat flow Qo 
goes to the cold source at To. 



The simplicity of the model allows the analytical eval- 
uation of other quantities Q . For instance, the energetics 
are straightforward. The mean heat flow released by the 
hot reservoir at T is 



(V/l + F)v, 



(5) 



since the system absorbs it to move against the external 
load F and potential force V/l. Therefore, the mean heat 
flow to the cold reservoir is 



(V/l)v, 



(6) 



because every time the potential barrier V is overcome, 
energy from the hot source is delivered into the bath 
at To. Finally, the mean power performed against the 
conservative force F is just 



W = Fv. 



(7) 



Note that the above expressions are consistent with the 
First Law, Q = Qo + W. Furthermore, the above charac- 
terization is compatible with the one introduced in Ref. 
0. See Fig. © for a scheme of the energetic quantities 
involved. 

According to the theory of nonequilibrium thermody- 
namics, the rate of entropy production can be expressed 



Q Qo 
T T Q ' 



(8) 



Now, this expression can be recasted in terms of a 
second order form involving the thermodynamic forces 
X x = F/T and X 2 = AT/T 2 , where the temperature 
difference is assumed to be small so that T — Tq = AT 
and T ~ To. Therefore, in the linear response regime we 
have 



S=(X U X 2 )[ ^ ^ 

\ -^21 ^22 



x 2 



(9) 



where Ly are the Onsager coefficients. We take Eq. 
© , substitute both expressions for the heats I© and © 
and expand the velocity v around the equilibrium state 
(F = and AT = 0) up to first order. This leads to 
an identification of the terms corresponding to every co- 
efficient in Eq. ©. After some tedious but standard 
calculations, one obtains 



L n =T 



sinh(a) 



VT ( a 
I \sinh(a) 



^22 



V 2 T 



sinh(cv) 



(10) 



(11) 



(12) 



where a = V/(2T). See |6j for a similar analysis. The 
Onsager coefficients provide a lot of information about 
the intrinsic nonequilibrium thermodynamic properties 
of the system. First one can check that the reciprocity 
relation, L12 = L21, is fulfilled and that the diagonal 
coefficients Ln and L22 are positive as they should. A 
relevant relation that is also found is 



L11T2 



T 2 

^12- 



(13) 



This makes the determinant of the L t j matrix equal to 
zero, which is in agreement with the Second Law. We 
can operate in a reversible regime and there is no en- 
tropy production. In fact, this condition is stronger than 
just telling us that, in the limit in which F = and 
AT = (motor at rest), no entropy is produced. More 
importantly, the above relation implies that the param- 
eter 



L12 



V L11L22 



-1. 



(14) 



which means that the motor is built with the condition 
of tight coupling This is the central result of this 

Report. 

When writing the efficiency of a thermal motor to low- 
est order in AT/T, it is found to be a function of the 
ratio of the thermodynamic forces k = X1/X2, so that, 




W 



AT 
~T~ 



_ K Ln\ K + L^/Ln _ (ig) 

L 2 1 J K + L22/ L21 



For the case in which |g| = 1, the last fraction of the 
above equation is equal to one. Now, if we stop both 
fluxes, J\ = J2 = 0, which would correspond to a van- 
ishing velocity and heat flow (but without the need of 
setting F = AT = 0), the above expression reduces re- 
duces to 



V 



1 



T ' 



(16) 
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which is the efficiency of Carnot. In general, for \q\ < 1 
it is always below Carnot's. The efficiency can also be 
found using the typical definition 



W 

? ? = "77 



Fv 



(17) 



Q (V/l + F)v V/(Fl) + V 
and then, by making v = 0, we obtain from Eq. J1J that 



T Fl 

n = 1 + V> 



(18) 



which inserted in Eq. I|17|) . gives again Eq. (|16[1 . For 

non-perfectly tight models, when stopping the power per- 
formed (at the stall force), there will still be a leak of heat 
which will make the efficiency far lower than Carnot's 
8J. It is true that achieving such upper bound in the 
efficiency of the motor in the limit of zero velocity has 
little practical relevance. In real devices it is more in- 
teresting to study efficiency at maximum power. It has 
been recently proved [9j that it is given by 



p/2 



(19) 



where p — q 2 /{2 ~ q 2 ). 
possible just when \q\ 
bound 



V 



C ]U! 

m 



Consequently, it is the highest 
1, leading to the Curzon-Alhborn 



1 



(20) 



The tight coupling of the present model ensures that this 
bound can be reached too. Let us check it by calcu- 
lating the efficiency in this regime. Using the relation 
amongst F, V , I, T and T derived from imposing that 
dW/dF = (the function found is hard to write explic- 
itly as a function of the force) , one can see that maximum 
power condition fulfills the condition 



V 

Yi 



(21) 



and thus, when inserted in Eq. (|17|) . leads to H20|) . which 
is precisely the Curzon-Alhborn efficiency for finite time 
endoreversible devices. In this case, one could perform 
finite time numerical simulations to recover such predic- 
tion, as it has been done for a similar model in Ref. Ijl. 



Finally, we want to address the question of converting 
the Brownian motor into a Brownian refrigerator |ll| . 
The system is acting as an engine when W > and 
Qo > 0. We define the refrigerator mode as Qo < and 
W < 0, since the effect of inserting work into the system 
leads to a heat flux leaving the cold source. Such oper- 
ation is sometimes regarded as a heat pump 12]. The 
existence of a refrigerator mode of operation is easily un- 
derstood recalling the previous linear irreversible thermo- 
dynamics analysis. To put it in words, since a thermal 
gradient drives the system against an external mechan- 
ical force, due to the cross coupling of the system {L 12 
and L21 coefficients), then an external mechanical force 
can lead to a thermal gradient opposing the one that al- 
ready exists. Therefore, any system that can work as an 
engine, has a region in the parameter space (which may 
be so small as to be hardy impossible to set) in which will 
work as a refrigerator. In Sakaguchi's model, as we can 
directly see, the reversal of the sign of the heat flowing 
to the cold reservoir happens simply when the velocity is 
inverted. Therefore, the transition from the engine mode 
of operation to the refrigerator occurs when the veloc- 
ity vanishes, and it is given by Eq. (|18fl . One technical 
comment to be added is that the fact that relation H13|) 
is found, leads to no separation between the engine and 
refrigerator areas 0- Obviously, the tight coupling also 
ensures an optimal performance of the refrigerator mode 
of operation. 

In conclusion, we have studied a simple model of ther- 
mal engine described by a Langevin equation which al- 
lows analytical calculations of relevant thermodynamic 
quantities. We give a fundamental explanation for the 
fact that the device can achieve optimal energetic perfor- 
mances. The question of whether a similar tight coupled 
construction is possible for a system in the underdamped 
regime is of crucial importance for the applicability to 
real microscopic machines, since heat flow driven by ki- 
netic e nerg y is relevant when inertia is taken into account 
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